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Abstract 

We introduce a method to construct conservation laws for a large class 
of linear partial differential equations. In contrast to the classical result 
of Noether, the conserved currents are generated by any symmetry of 
the operator, including those of the non-Lie type. An explicit example is 
made of the Dirac equation were we use our construction to find a class of 
conservation laws associated with a 64 dimensional Lie algebra of discrete 
symmetries that includes CPT. 



1 Introduction 

The connection between symmetry and conservation laws has been inherent 
in all of mathematical physics since Emmy Noether published, in 1918, her 
hugely influential work linking the two |20| . Since then many have put forward 
approaches to study conservation laws, through a variety of different means 
[2 O [TOl [Til [15]. In each case, a conservation law is defined as follows. 

Definition 1. Let A[u] = be a system of equations depending on the indepen- 
dent variables x — [x^ , . . . , a;"), the dependent variables u — {u^, . . . , u™) and 
derivatives thereof. Then a conservation law for A is defined by some P — P[u] 
such that: 

DivF = (1.1) 

A=0 

where [u] denotes the coordinates on the iV-th jet of u, with N arbitrary. 
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Noether's theorem is applicable in the case where A[u] = arises as the 
Euler-Lagrange equation to an associated variational problem. It is well known 
[HI [211 131 13 that a PDE has a variational formulation if and only if it has self- 
adjoint Prechet derivative. That is to say: if the system of equations A[u] = 
is such that Da — DJ^ then the following result is applicable. 

Theorem (Noether). For a non-degenerate variational problem with L[u\ = 
J^£,dx, the correspondence between nontrivial equivalence classes of variational 
symmetries of L[u\ and nontrivial equivalence classes of conservation laws is 
one-to-one. 

Here the symmetries we speak of are generalised symmetries which corre- 
spond to those with generators of the form: 



Without loss of generality, one may confine attention to generalised vector fields 
in evolutionary form, where: 



The theorem of Noether then provides equivalence classes of conservation laws 
in terms of the characteristics of Xq, i.e the Qi[u]. It can be shown '21\ that 
the Qi are related to the corresponding conservation laws via: 



which holds for all u. This clearly gives rise to a conservation law since the LHS 
vanishes on solutions. Here we say Q is the corresponding characteristic to the 
conservation law. 

More generally, for systems that do not obey a variational formulation, we 
can consider (|1.2p the definition of a conservation law. It can be shown [21] that 
the condition is equivalent to the existence of some Q[u] such that Q • A is 
a total divergence, as in ()1.2|1 . Since the Euler operator annihilates divergences, 
this condition reads: 





Q • A = DivP 



(1.2) 



E,{Q-A)=0 l<i<m 
for each {x,u). Here E denotes the Euler operator: 



(1.3) 
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where a is a multi-index and it^ = We regard (|1.3p as the defining 

property for Q to be a characteristic for a conservation law for A. In this case, 
there is no interpretation of the Q in terms of symmetries of the equations. 
Indeed, expanding the condition in (|1.3p gives: 

Da(Q) + D^(A)-0 

for each (x, u). In particular, on solutions we have what Anco and Bluman [I1I2] 
refer to as the adjoint- symmetry condition: 

D^(Q) =0. (1.4) 

A=0 

This gives us a necessary, but not sufficient condition for Q to be a characteristic 
for a conservation law. Now if A[u] =0 corresponds to a variational problem, 
then — Da and (|1.4p is automatically satisfied if Xq is the infinitesimal gen- 
erator of a symmetry of A[m] = 0. This has an obvious connection to Noether's 
theorem: if Xq corresponds to variational symmetry of a variational problem, 
then it is also corresponds to a symmetry to the corresponding Euler-Lagrange 
equations (the converse is not true in general). 

Throughout our discussion regarding symmetries, it has been implied that 
these symmetries correspond to a locally connected Lie group of transforma- 
tions. If we were to talk of more general types of symmetry, then much of the 
previous material becomes void because of our use of infinitesimal generators. 
In addition, in our general definition of a conservation law, we insisted that it 
be a total divergence depending only on [u] - if we relax this, then it allows for a 
larger class of nontrivial total-divergences that vanish on the space of solutions. 
In this paper we address these points for a class of linear PDEs. 

We consider a general set of symmetries, referred to as non-Lie symmetries 
[19l [iHl [13] . This class includes generalized symmetries, non-local symmetries 
and discrete symmetries. An exposition into the construction of non-Lie sym- 
metries can be found in [TH] where several examples are given. Given that this 
class of symmetries is far larger than those considered in the classical work of 
Noether, there is potentially an even stronger correspondence between symme- 
try and conservation laws for PDEs. 

As an example, consider the scalar field u{t, x), x dVl C R" whose evolution 
is prescribed via the heat equation, ut = Au. Then the following is a nontrivial 
conservation law for this system: 

Dt [u^'u] + Div [u'^Vu - uVu"] = 
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where u*(a;, t) = u{x, s — t), s fixed. Taking, for example, the solution to the re- 
lated Cauchy problem with initial data / e gives the family of conserved 
quantities 

= [47rt{s-t)]-/^ ^^^e-ll--'^ll^/-^*-|l---ll^/^(--*)/(y)/(^)d"xd"yd"^ 

for t e (0, s). This result can not bo achieved using standard techniques. 

We formulate a method that provides insight into constructing new con- 
servation laws generated from non-Lie symmetries, and that does not require 
a variational formulation. The method returns the classical results when the 
symmetries are of the Lie type, but provides new, non-trivial conservation laws 
when algebras of discrete symmetries are used. We highlight this result by con- 
structing new conservation laws for the Dirac equation, corresponding to a large 
Lie algebra of discrete symmetries. For example, we show that CPT symmetry 
in the Dirac equation generates the operator (after quantisation): 

which can be shown to commute with the Hamiltonian of the theory. 



2 Construction of the conservation law 

We work on a formal level so that all functions in question have sufficient 
smoothness and decay so that relevant integrals are well-defined. In particu- 
lar, for evolution equations on (0, T) x C R"+^ with conservation laws of the 
form: 

^ + Div J = 
at 

where p = p[u] and J = J[u], we assume the solutions are such that J^pdx con- 
verges and J\dQ = 0. Throughout we work with functions of n + 1 independent 
variables so C°°(R"+\R'") etc. wiU be denoted simply by C°°(R™). We con- 
cern ourselves with linear PDEs, so that in our previous notation A[u] = C[u] 
where £ : C°°(R'") C°°(R') is a linear matrix differential operator with 
smooth coefficients. We denote the formal adjoint of £ hy C* : C°°(R') — > 
C°°(R'"), the coefficients of which are defined by (£*),;; = (Cji)*, meaning that 
one first takes the transpose of the matrix differential operator £, then takes 
the formal adjoint of each of differential operators. 
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When we talk of symmetries of a PDE, we refer to them in the sense of 
Fushchich and Nikitin 14]. The definition we use is as foUows. 

Definition 2. We say the operator F is a symmetry of the linear PDE A[u] = 
C[u\ =0 if there exists an operator ar such that: 

[C, F] - arC 

where [•, ■] denotes the commutator by composition of operators so CT = CoT. 
We denote the set of all such symmetries by sym(A). 

There is a rich source of literature on the study of such symmetries of linear 
PDOs [13l [Ml [TBI dl] some of which we shall use in the sequel to construct 
conservation laws for the Dirac equation on Minkowski space. The analysis of 
such symmetries is not confined to those mentioned earlier, but also to those of 
the non-Lie type including discrete symmetries. 

Remark. Definition [2] can be viewed in terms of recursion operators. Recall 
that 7?. is a recursion operator for the PDE A[u] = if there exists an TZ such 
that: 

D^n = 7^DA. 

In our case, A[u] is linear, so Da coincides with C In terms of the notation in 
definition [5] we have T = TZ and ar = TZ — TZ. 

There is also a correspondence between these symmetries and the usual Lie- 
type as found in the standard treatment. Indeed, suppose the linear, A'^-th order 
PDE is invariant under the action of a one parameter group of transformations 
Fe = exp(eX), with corresponding infinitesimal generator: 

i—1 i—1 



Then by definition of X we must have: 

pr^(X)(AM) 



= 

A=0 



where pr^ (X) is the prolongation of X onto iV-th jet space. For this to hold, 
it must be the case that (see [21] p. 82): 

M 

pr^(X)(AM) = ^^^,(x)AAM, 

1=0 

for some Fi{x) that are unimportant and some M > 0. Now the LHS is equal 
to Da{Xu), and since A is linear, this is just C[Xu]. From this we see that the 
relevant commutator [C,pr^ (X)] takes the required form. 
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From here on we deal entirely with symmetries of the form in definition [5) 
An immediate consequence of this definition is as follows. 



Lemma 1. For each T £ sym(£) and u G ker(>C) we have Tu £ ker(Z!). 

We saw in fJT] that the existence of a conservation law to the equation 
A[ii] =0 was dependent on finding solutions to the so called adjoint-symmetry 
condition: 

Dl(Q) - 0. (2.1) 

A=0 

In the linear case with A[u] = C[u\, this condition reads: 



C* 

The following lemma makes the condition also sufficient 



= (2.2) 

A=0 



Lemma 2. For a linear PDF A[u] — C[u\ ~ Q, there exists a conservation law 
for each solution to 



Proof. We present an argument similar to that presented in [23j. Given £, we 
associate with it the bihnear map He : C°°(R™) x C°°(R') C°°(R), where 
is defined explicitly by: 

Ilc{Q,P)^Q-C[P]-P-C*[Q] (2.3) 

where P, Q are m-tuples (respectively ^-tuples) of smooth functions. Integration 
by parts shows this is a total divergence for each P, Q, so that Iic{P, Q) = DivX 
for some X depending bilinear ly on P and Q. Alternatively one can view this 
as a consequence of Green's formula [23]. Setting P — u we see that: 



DivX|A=o = -u-C*[Q] 



A=0 



If Q solves (I2.2p then DivX = on solutions to A[u] — 0, i.e a conservation 
law. ■ 

Corollary 1. If C is self-adjoint or skew-adjoint, then each T G sym(£) gen- 
erates a conservation law. 

The corollary follows by setting Q = Tu and appealing to the result in lemma 

m 
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Example. Consider the wave equation on (0, T) x R", £ = □. Constructing 
the bihnear map as in the proof to lemma [2] gives: 

UciQ, P) = DtiQPt - PQt) + Div (FVQ - QVP) 

Setting P — u and Q = Tu with T G sym(n) gives a conservation law. Integrat- 
ing over R" and the divergence theorem then gives the conserved quantities: 



I [(ru)uf — uDtiTu)\ Ax 



for each symmetry F of □, assuming that all fields and their derivatives have 
sufBcient decay at spatial oo. Choosing T = dt gives rise to conservation of 
energy: 



R" 



(uj - uutt) dx 
fu^ + lVuP) dx, 



which follows from integration by parts. This is in agreement with the classical 
result. 

Example. We consider the coupled, linear KdV-KdV system [S1[T2] on (0, T) x 
R: 



Ut + Vxxx + — 0, 
Vt + Uxxx + Ua; = 0. 

Writing the equations, A[u,w] = 0, in matrix form gives: 



(2.4) 
(2.5) 



Dt Dl + Dx 
Dl + Dx Dt 





u 




V 



= 



We note that the system is skew-adjoint, and so the preceding results give a 
conservation law for each F G sym(A). Integrating this conservation law over 
R gives the conserved quantities: 



Kr 



R 



[7i(Fu) 4- w(Fu)] dx 



for each F £ sym(A). Computing the symmetries by hand is straight forward, 
using the methods of Lie. The following vector fields correspond to the genera- 
tors of Lie-point symmetries with linear coefficients: 

Vi = dx, V2 = dt, V3 = du, V4 = dy, 

V5 = udv + vdu, Ve = vdu + udv, V7 = tdv - xdu, Vs = xdv ~ tdu 
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which generate an 8-dimensional Lie algebra. Labelling the corresponding con- 
servation laws Ki we find ki, K2 are trivial and K3, K4 follow from the fact that 
(|2.4p - (j2.5p are already in conservation law form. Then K5, Kg corresponds to the 
conservation laws: 



K5 = / {u +v )dx, Kg— uvdx 
Jr Jr 

Finally, the vector fields Vj, Vs generate the conserved quantities: 

kt — {tu — xv) dx, Ks = {xu — tv) dx 
JR Jr 

Now we make use of a family of discrete symmetries defined by Ts(u,v) = 
{—u", ti^) with v''{x, t) = v{x, s — t) etc, with s fixed. Then for t G (0, s) we have 
the conserved quantity: 

Kg — [vv'^ — uu^] dx. 
Jr 

To explicitly show this is indeed conserved, we compute the derivative. 



dK f 

"d^ " /r ""''"'^''"^ '■"^ Kxx)v + {vxxx + v^)u'' - {v^^^ + vl)u] dx 

Dx [vu'^ — MW' + u'^Vxx + Vu'^xx ~ "x^a: ~ V^Uxx — UV^^.j, + U^V^] dx 

R 

= 

which holds for each t E (0, s). 

In this example we used an obvious discrete symmetry of the system. In 
general, to find the discrete symmetries of a system of PDEs we can apply the 
result of Hydon [16j . In this method, one considers a general symmetry: 

F : {x{x,u),u{x,u)) 

that may be continuous or discrete, then "dresses" this symmetry with a one- 
parameter group of symmetries Fi(e) = exp(eXi). The resulting symmetry: 

f,(e) = FF,(e)F-i ^ exp(eFX,F-i) 

is of the one parameter form, generated by Xi — rX,F~^. The determining 
equations for Xi can then be related to the Xi, and the result is a system 
of nonlinear constraints that can be solved to give the most general form of 
symmetry. In particular, the discrete symmetries can be found by quotient- 
ing out the symmetries found via Lie's method. Several worked examples are 
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given in [16] which demonstrate the algorithmic process involved: the number 
of computations is not much more than that needed to compute the continuous 
symmetries. 

We now turn our attention to the study of systems of PDEs with constant 
coefficients, and provide a general result that generates conserved quantities for 
each symmetry of the system. First we define the type of linear operators in 
question. 

Definition 3. We denote the space of linear operators spanned by M"_Dq,, 
Af„ e Mat,„(C), by L^(C"). If £ £ L=(C™) then £ : C°°(C'") ^ C°°(C™). 



Such PDEs include both the wave equation and the linearised KdV-KdV 
system discussed earlier. Another example would be the linearised Navier-Stokes 
equations [1]: 



Ut + Vp — v/S^u = 
Div(u) = 



(2.6) 
(2.7) 



where u = u{x, y, z, t) is the velocity field, p = p{x, y, z, t) is pressure field and 
V > {). Then the corresponding element of is given by: 



L = 



Dt -i^A 
Dt-i^A 













Dy 

Dt-vA D, 




and equations (|2.6p - (|2.7p are given by £[(/)] =0 with 4> = {u^, u^,u^,p)- 

Definition 4. We say £ MatTO(C) are simultaneously semi-conjugate to 
[M"]t if 3Ai,A2 e GL{m,C) independent of a such that [M°]t = AaAf^Aj-^ 
for each a. 

Lemma 3. Suppose C = M^Da £ L'^{C"^) and the M" are simultaneously 
semi- conjugate to [M"]^ . Then then there exists a pair of invertable linear 
tors, Pi,P2, such that C* = 'P2CV^^ . 



Proof. Using C* = Ea(~l)'"' find: 



(2.; 
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Now we introduce the parity operator P : C°°(C'") C°°(C™) defined by 
P(j>{x) — x). We make the following observation: 



Using this identity in (|2.8p and the linearity of the associated operators we find: 



Then from the definition of C we see the required operators are 7^2 = A2P and 



Here we see the parity operator playing an important role, and it needs 
clarification as to whether this will cause any problems with our analysis. Terms 
such as P[(j)] will be present in our conservation laws, and as such care must be 
taken when forming the corresponding conserved quantities. 

In our opening remarks, we have made the assumption that in the case of 
evolution equations on (0, T) x J7 C R"+^ all fields are such that the relevant 
"flux", J, vanishes on dfl. In most cases this requires that the fields have 
sufRcient decay at spatial 00. If we say, for example, that •) g W'''P{fl), 
then it is clear that the discrete symmetries on the spatial independent variables 
will not cause any problems, assuming they are well-defined. For example if 
n = R", then if u{x,-) G W'''P{n) then u{-x,-) £ W'^'Pi^) also. Similarly if 
is any isotropic domain. However, in the case of discrete symmetries of the 
temporal variables, they must take the form F : u{-,t) u{-,s — t), valid for 
t G (0, s) and s G (0,T). This means we need to modify our Parity operator so 
that when it acts on temporal variables, it does so in this prescribed form. As 
we shall see, however, in many instances the parity operator can be dispensed 
with altogether, by modifying the matrices Ai, A2 in (|2.9p - particularly in the 
case of Dirac operators (see lemma [5] and corollary [2]) . 

In what follows we provide some lemmas that provide sufficient conditions 
for the matrices M" to be simultaneously semi-conjugate to the [Af"]'''. 

Lemma 4. Suppose L = J^a^^^D^ G L=(R™) such that [Af",M'3] = for 
each a, p. Then C* = V2'0.Pi^ for some linear, invertable T'i,P2- 



= P[D^{P^)] 

^{PoDaOl 




(2.9) 



Pi = p^-i = A-^P. 



10 



Proof. Since the M" commute we may choose a basis {e^}™ such that the 
M" simultaneously assume Jordan canonical form. This gives — SiJ^S^^ 
where J" is in canonical form and 5*1 G GL(to,R). Now introduce a new basis 
such that: 

Bi Bi ^ ejn-i+i, 1 < i < m 

Under this change of basis, the J" [J"]*. This then gives M" = S'2[J"]*52"^ 
for some ^2 £ GL{m, R). The result follows since [J"]* = 6*3 [A/"]* 6*3"^ for some 
5*3 e GL(to,R). Then £* = V2CV^^ for some linear, invertable Vi,V2- ■ 

Lemma 5. Suppose C — A/^Dq, G L'^iTU"') such that the Ma belong to a 
unitary representation of the Clifford algebra C^i^g(R). Then C* ~ T'2jO.Vi^ for 
some linear, invertable 'Pi , 7^2 ■ 

Proof. Recall that if {7'} are a unitary representation of C^i^g(R) then they 
necessarily obey the anti-commutation relation: 



From this and the fact that 7' [7*]^ = [7*] ^7* = I we find: 

[-7', l<i<q + l 

From this we deduce the result: 

[7']t^-^iy^\ l<i<q+l. 

So each of the M" are simultaneously semi-conjugate to the [M"]'''. The result 
now follows from lemma [3) ■ 



We now provide a concrete example to illustrate lemmas [3] and HI 
Example. Consider C — J2a M"Da G L^(R^) given by: 

13 



£.= 



Dt + D'l DtD, 
Dt + Dl 





1 







1 










1 






A + 













1 







1 











Clearly each of the M" commute, and since each matrix is already in Jordan 
canonical form we only need to compute the appropriate matrix S2 G GL{2, R) 
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as in the proof to lemma [H Using the result in (|2.9p gives: 



£* = 



"o l' 




"i o' 




"i o' 




"o l' 




l' 




'( 




Dt + 








DtD^j P 




1 




1 




1 







1 



1 


o' 




1 


o' 


Dl + 





o" 






Dt- 













1 







1 




1 






DtD, 



-Dt-Dl 
DtD, -Dt~Dl_ 

Which agrees with the expected result. 

We note that there exists an alternate argument for the result in Lemma [3] 
which applies to all elements in L'^(R"'). Indeed, if we consider £ as a matrix 
over the ring of pseudo-differential operators O [22] of the form: 

finite 

X.v=Y^ c"'{<l>)Da (2.10) 

— oo 

we have considerable freedom. In our case, the c" are constants and as such 
the associated pseudo-differential operators form a commutative ring. What is 
more, inverses exist - so on a formal level we may treat our C as matrices over 
this field of differential operators. Then recalling from linear algebra that each 
matrix over a field F is similar to its transpose, we may use a similar argument 
to that found in lemma[3]to prove the result for arbitrary C € L'^(R). However, 
as we shall see, we will be interested on the action of the relevant operators 
Vi on test functions <j) G C°°(R'"), so the inclusion of operators of the form in 
(|2.10p causes considerable technical issues of how one should define [(p] . For 
this reason we treat the indicated subset of L'^(R™). 

Theorem 1. Let A[u] = &e a linear system of PDEs associated with the 
operator C — M" Da G L'^(R™), where the M" are simultaneously semi- 
conjugate to the [Af"]^. Then for each T G sym(A) there exists a conservation 
law. 

Corollary 2. Suppose the C associated with A[u] = is a polynomial in the 
Dirac operator p = ^.^YDi. Then for each T G sym(A) there exists a conser- 
vation law. 

Proof. We see from lemma [1] that for each F G sym(A) we have Tu G ker(£) 
for solutions to A[ii] = 0. Also by lemma [3] it follows that if w G ker(£) then 
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Vi G ker(£*). Indeed, using the result from lemma [3] we have: 
CV^^u = V2CViV^^u = V2CU = 

since u G ker(i2). Combining these results gives: 

u G ker(£) ^ p-^Tu G ker(i;*) 

for each F G sym(A). This can be restated as: 



= 

A=0 



for each F G sym(A). The result now follows after an application of lemma [2j 
The corollary follows from this result and an application of lemmas [3] and [3 ■ 

We give an extended example in the next section which utilizes the results 
in theorem [1] and corollary [21 



3 The Dirac Equation 

Dirac's celebrated equation governing spin i fermions takes the covariant form: 

[ip-m)^ = Q, p = -l'^Dt-"f'D^ (3.1) 

where here and throughout this section we employ summation convention (ro- 
man letters running over 1,2,3 and Latin over 0,1,2,3). The {7^} are members 
of the Dirac algebra Cli,^{C) = C/i^3(R) ® C which obey: 

{^f , 7^7" + 7"7^ = Syy^'-I (3.2) 

where 77 is the metric on Minkowski space, with signature (+, — , — , — ). Here '0 
is a Dirac spinor associated with the (^, 0) © (0, \) representation of the Lorentz 
group. We employ the Dirac representation of C^i^3(C), so that the {7^} are 
given by the following 4x4 matrices: 

where the cr' are the Pauli spin matrices. We use r]^'^ to raise and lower indices 
so that 7^ = Vfii^l'^ where rj^,j'q'^P = 5'^. The field equation (|3.ip comes from 
varying the following action: 

= / — ni)^p dx. 
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£ = 



(3.4) 



Here = ^p'''^^ is referred to as the Dirac conjugate. This variational formula- 
tion admits a J7(l) variational symmetrjQ, from which we deduce the conserved 
current = tpj^tlj via Noether's theorem. We will now derive this result with- 
out referring to variational principles. We identify the appropriate £ G L^(R'*) 
as: 

(iDt - m)l2 -iaWi 
ia'Di -{iDt + m)l2 

Now since this is a (linear) polynomial of Dirac operators, corollary O applies 
and we may construct a conservation law from each symmetry of p.ip . Noting 
that the fields are complex, we construct n£ using the adjoint associated with 
the inner product (/, g) — J fg* dx. A routine calculation gives: 

Since £ is a polynomial in ^, the result in lemma [5] gives: 

£* =7"(^^-m)7° (3.5) 

Then in the notation of lemma [Sj we have that Vi ~ 'P2 =7"- Now using the 
result of theorem [T] with F = id G sym( A) gives the conservation law: 

n£(v,^) = 5^(*(7V)V^) 

agreeing with the classical result corresponding to conservation of charge. How- 
ever, we may now use the result from the previous section to create more con- 
servation laws. If r G sym(A) then by theorem [TJ 

H£(^,7°rV') = (3.6) 

is a conservation law associated with the Dirac equation. Each of these generate 
a conserved quantity: 

Kr = y (r?A)Vda; (3.7) 

which follows from l|3.6p and (7")^ — I. This holds for all symmetries of (|3.ip . 
not only those corresponding to variational symmetries in the Lagrangian for- 
mulation. For completeness, we observe that the result in (|3.7I) agrees with the 
classical results when F G sym(A) is a symmetry of the Lie type. 



^Or gauge symmetry in the physics literature. 
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It is well known that the maximal algebra of Lie type symmetries of the 
Dirac operator is the Poincare algebra, which has basis: 

We focus on the symmetries that generate rotations, i.e: 

It is a simple exercise to show that this gives rise, by means of Noether's theorem, 
to conservation of angular momentum (orbital and spin) J where: 

J = y"v^t (xAp+ iS)i/'dx (3.8) 

where p = — iV and S = diag(tTi, tXi). We use the same subset of the Poincare 
algebra that generates rotations, so the associated symmetry operator is: 

Ti = ieyfc {xjPk - XkPj + 
= e,,k^,pk + 

We then see this leads to the same conserved quantity as derived in p.8p . Clearly 
each of the conserved quantities that follow from Noether's theorem can also 
be generated by this method, but the important point to outline is that we 
may construct other conserved quantities, generated by non-Lie symmetries, 
including discrete. 

In [18j the algebraic aspects of the discrete symmetries of the Dirac equation 
were studied in detail. We follow this work, and study the algebra associated 
with the discrete parity, charge and time symmetries present in the Dirac equa- 
tion. The basis for this algebra is: 

= = 1746*, Fs = i"f2C, Tq = iPs 

where 02'>P{^) = V-'C^^Oj^;!, — X2,a;3) etc, 9ip{x) = ip(—x), dip = ip* ^-nd 74 is the 
fourtli gamma matrix: 

74 = 170717273- 

These basis elements actually correspond to a 7-dimensional Clifford algebra: 
{Pa, Pb} = gabl, 9 = diag(+l, -1, -1, -1, +1, +1, +1). 

^Or often referred to as the fifth. 
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It can be shown that the enveloping algebra of this ClifFord algebra is isomorphic 
to 0[(8, R) which means there are 64 linearly independent combinations of these 
elements. The algebraic aspects of these symmetries has a rich structure, but 
we shall not pursue it here so refer the reader to jT31 HU [TS] . We now form 
a conserved quantity from this an element of this algebra: we choose T = 
Tq G sym(A) for simplicity. It follows from (|3.7p that the following quantity is 
conserved: 

= ■ip{-t,x)j4^p{t,x)dx. 



dt 



Noting that ip{—t,x) satisfies idoipj'^ = —idit/jj'^ + niip, we find: 

[idiip{-t,x)j'j^'il!{t,x.) - i'0(-t,x)7'*7'9,7/>(t,x)] dx 
[td,^Pi-t, ^h'-f^^P{t, x) + ^^(-^, ^)fj%,p{t, x)] dx 
[V;(-i,x)7'7V(i,x)] dx 



= 0. 

Where we used the fact that 74 anti-commutes with the 7^ and the assumption 
that the fields decay at spatial 00. We can express this quantity in terms of the 
usual creation and annihilation operators - a short calculation reveals: 

'^o - / E {--pb^ + b%'a;) (3.9) 

where we have employed the usual quantisation procedure of promoting the 
fields to operators on a Fock space, and the a*, etc. obey the anti-commutation 
relations: 

in the distributional sense, and all other anti-commutators being zero. We see 
Ko annihilates the vacuum, and using the usual anti-commutation relations we 
find: 

So kq swaps an anti-particle for a particle, and reverses the direction of mo- 
mentum. A routine calculation shows that kq commutes with the Hamiltonian 
H: 

[H, «o] = 0, H = J^E^Y1 + 

so it is indeed conserved. It is interesting to note the conserved quantity associ- 
ated with the CPT symmetry that is inherent in all field theories governed by a 
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Hermitian Hamiltoiiian. In this case it is given by r4r5, and the corresponding 
conserved quantity is: 

K45= yd3a;i^(-i,-x)*7V7V(i,x), (3.10) 

which is easily shown to be constant in time. Expressing this in terms of the 
creation and annihilation operators for the particles and anti-particles, we find 
that the operator in (|3.10p is equivalent to: 

/(§^E(-ir(«pV^ + ^;^^'^0' (3.11) 

where the spin sum is over Z/2Z. The derivation of this result can be found in 
the appendix. A routine calculation confirms that this operator commutes with 
the Hamiltonian. 



4 Conclusions 

In this paper we have seen that the number of conservation laws for a system of 
linear PDEs can be greatly increased if one considers the associated densities to 
depend not only on [u], but also discrete transformations of thereof. Allowing 
for this extra freedom, we have produced a general result applicable to a large 
class of linear PDEs, including those associated with polynomials of Dirac op- 
erators. The results of Hydon [16] and Fushchich [13] allows for an algorithmic 
construction of the non-Lie symmetries, upon finding which theorem [T] provides 
us with a means to construct corresponding conservation laws. 

Acknowledgments. The author would like to thank Thanasis Fokas for his 
lessons in perseverance and George Weatherill for many useful discussions. In 
addition, the author thanks the anonymous referee who made him aware of the 
work in This work was supported EPSRC. 
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5 Appendix 



CPT Symmetry in the Dirac Equation 



It was shown in §3 that the CPT symmetry inherent in the Dirac equation gives 
rise to the conserved quantity: 



(5.1) 



where the spin sum is over Z/2Z. In this section we outUnc the details of the 
calculation that produces this result, after applying the standard quantisation 
procedure to the term: 



d^x V(-i, ^Kt, x), 



(5.2) 



The calculations in this section would be somewhat easier if we chose to work 
with Marjona spinors from the offset, but we shall stick with our initial choice 
and suffer the consequences. Firstly we make an observation: 



-i [7*]^ 9p — m)tlj{xy 



+ij"dt - ij^dx + ij^dy - ij^d^ ~ m)i/'(x)*, 
ip - m)il;{-t,-x,y,-z)* 





^ = 
= 
= 
= 

And so we set: 

'4;{-t,-x)* = i}{t,x,-y,z). (5.3) 

We shall use this in (|5.2[) . along with the mode expansions for i/j, -ip. In the 
Schodinger picture the relevant mode expansions are given by: 

where p' — {p^ , —p'^ , p^) . Using these expressions in (|5.2p and doing the x 
integral and one of the p integrals, we find: 

d^p 1 



(5.4) 
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Using the anti-commutation relation for the gamma matrices in (|3.2p . we have 
^2^o^4 _ ^o^4^2 jnggj-ting the 7^ matrix into the right hand side bracket in 
(|5.4p leaves us to compute j^Us{p') and 7^Ws(— p')- We consider 7^Ms(p') only, 
since the computation for 7^Ws(— p') is entirely analogous. In our representation, 
7^1*3 (p') is given by: 



a' 

-a^ 



. E„,+m 



where xi — [J] ^^'^ X2 — [?] are 2-spinors. Using the fact that Epi — Ep and 
tT^(CT ■ p') = — (cr • p)cr^ we find: 



/ (°-p)xA 

727.,(p') = *(-l)^t>,+i(p), where i;,(p) = V^p' + ( 1 

where s £ Z/2Z. A similar computation for 7^Ws(— p') reveals: 

7V(-p')=*(-l)^+^.+i(p). 
The expression in (|5.4|) now takes the form: 

X 7V (a;,(-l)^i;.+i(p) + 5ip,t(_i)s+i^^^^(p)^ , (5,5) 

Now we must consider terms of the form 'itr(p)7'^7*fs+i(p) and similar. We 
perform the calculation for Ur{p)'^'^'y^Vs+i{p) since the computation is similar 
for the remaining three terms. In our representation the computation reads: 



Ur[P 



(p) ( X ) Vs+i{p) = (Ep + m) (xj, -xl 



\ i'^-p) 1 


)( 


Ep+m 





Xs+1 

('^■p) 



{Ep + m) 

[Ep + m + 
2EpS'"''+^ 



xlxs+i + xl 



-Ep+m 



E„+m 



Xs+1 



Xs+1, 



Ep — to] S 



■r.s + l 



pp = Ep — m^. A similar calculation reveals: 



in which we used (cr • p)^ 

Vr{-Ph°l^Us+i{p) ^ 2EpS^''+\ 

with the remaining two terms in (jS.Sp evaluating to zero. Using these results 
we find that (|5.5p reduces to: 

(5.6) 



/(|;^E^(-ir(«;V^-^V-"p' 
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Now discarding the unimportant constant z, changing p > — p in the second 
term and using the anti-commutivity of the &'s, we recover the expression in 
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